Gluon shadowing and hadron production in heavy-ion collisions at LHC 
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The recently published first measurement of charged hadron multiplicity density at mid-rapidity 
dN ch /dri = 1584 ± 4(stat.) ± 76(sys.) in central Pb + Pb collisions at y/s = 2.76 TeV by the ALICE 
Experiment at LHC is in good agreement with the HIJING2.0 prediction within the experimental 
errors and theoretical uncertainties. The new data point is used to carry out a combined fit together 
with the RHIC data to reduce the uncertainty in the gluon shadowing parameter s g which controls 
the overall magnitude of gluon shadowing at small fractional momentum x in HIJING2.0 model. 
Predictions on the centrality dependence of charged hadron multiplicity density at mid-rapidity with 
reduced uncertainties are given for Pb + Pb collisions at y/s = 2.76 and 5.5 TeV. The centrality 
dependence is surprisingly independent of the colliding energy similar to that in Au + Au collisions 
at RHIC for most of centralities starting at Af part = 50 (100) at y/s = 2.76 (7) TeV. However, it 
becomes stronger in peripheral collisions at higher colliding energies. 
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Bulk observables such as rapidity density of hadron 
multiplicity and transverse energy provide important in- 
formation on the initial entropy and energy production 
in high-energy heavy-ion collisions. They also provide 
constraints on the initial conditions for hydrodynamical 
study of the collective phenomenon and hard probes such 
as jet propagation and suppression. Because of the non- 
perturbative and many-body nature of the physics pro- 
cesses involved, a first principle calculation of the bulk 
hadron production is so far unaccessiblc in Quantum 
Chromodynamics (QCD). Instead, one has to reply on 
theoretical and phenomcnological models to estimate the 
bulk hadron production in high-energy heavy-ion colli- 
sions. Even with constraints by experimental data from 
the Relativistic Heavy-ion Colliders (RHIC) experiments 
[l[ there still exits large uncertainties in the theoreti- 
cal and phenomcnological estimates of charged hadron 
multiplicities in heavy-ion collisions at the Large Hadron 
Collider (LHC) both among different models and within 
each model Q . 

Recently, ALICE Experiment at LHC published the 
first experimental data on the charged hadron multi- 
plicity density at mid-rapidity in central Pb + Pb col- 
lisions at y/s = 2.76 TeV 0). The measured dN ch /drj = 
1584 ± 4(stat.) ± 76(sys.) for the top 5% central Pb + Pb 
collisions is larger than most of theoretical and phc- 
nomenological predictions, including all the latest color- 
glass-condensate model calculations Q. Such an unex- 
pected large hadron multiplicity will have important im- 
plications on the underlying mechanism for initial parton 
production. It will also have important consequences on 
the study of other phenomena such as collective flow and 
jet quenching in Pb + Pb collisions at the LHC energies 
since they all depend on the initial condition for the bulk 
matter evolution. 

The first ALICE data pi is in good agreement with 
the HIJING2.0 prediction |1 @ within the experimen- 



tal errors and theoretical uncertainty which is controlled 
mainly by the uncertainty in the parameterization of the 
unknown nuclear shadowing of gluon distribution. The 
new HIJING parameterization of the gluon shadowing in 
nuclei Q was constrained mainly by experimental data 
on charged hadron multiplicity and its energy and cen- 
trality dependence in heavy-ion collisions at RHIC within 
HIJING2.0 Q model which is an updated version of the 
original HIJING1.0 model @, Q. At the LHC energies, 
initial parton production probes gluon distribution at 
much smaller fractional momentum x. The range of the 
gluon shadowing parameter s g = 0.17 — 0.22 allowed by 
the RHIC data, which controls the overall magnitude 
of the gluon shadowing at small x in the new HIJING 
parametrization, leads to much larger uncertainties in 
the final charged hadron multiplicity, up to about 15% 
in the most central Pb + Pb collisions at y/s = 2.76 TeV 
@. Within the same HIJING2.0 model, the ALICE data 
provides a much stringent constraint on the gluon shad- 
owing. In this note, we will carry out a global fit of 
the combined RHIC data on the centrality dependence 
of charged hadron multiplicity in Au + Au collisions and 
the new ALICE data in the most central Pb + Pb colli- 
sions at LHC to provide a better constraint on the gluon 
shadowing parameter s g . With a smaller range of the 
gluon shadowing parameter s g = 0.20 — 0.23, we will pre- 
dict with reduced uncertainty the centrality dependence 
of charged hadron multiplicity density in mid-rapidity for 
Pb + Pb collisions at both y/s = 2.76 and 5.5 TeV. 

HIJING [f| 0] is essentially a two-component model 
for hadron production in high energy nucleon [8hJl0j and 
nuclear collisions [U, [T^ • In this two-component model, 
one divides nucleon interaction into soft and hard part 
separated by a cut-off po in the transverse momentum 
transfer between colliding partons. Jet production with 
transverse momentum pt > po can be calculated within 
the collinear factorized perturbative QCD (pQCD) par- 
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ton model while the soft interaction is characterized by a 
parameter in the effective cross section a so f t . These two 
parameters, <J S oft and po in HIJING, are determined phc- 
nomenologically by fitting the experimental data on total 
cross sections and hadron multiplicity in p+p/p collisions 
[lijj . HIJING2.0 [H is an updated version in which old 
Duke-Own (DO) parameterization [l4| of parton distri- 
bution functions (PDF's) is re plac ed by the Gluck-Rcya- 
Vogt(GRV) parameterization [15(. Because of the much 
larger gluon distribution in GRV than DO parameteriza- 
tion at small x, one has to assume an energy-dependent 
cut-off po(\fs) and soft cross section a so jt{\fs) @ i n order 
to fit the experimental values of the total and inelastic 
cross sections oip+p/(p) collisions. The values of po and 
a so ft are further constrained by the energy-dependence 
of the central rapidity density of the charged hadron mul- 
tiplicities. This updated version HIJING2.0 can describe 
most of the features of hadron production in pp collisions 
at colliding energies up to 7 TeV at LHC 131. 

For high-energy heavy-ion collisions, both nuclear 
modification of the parton distribution functions and jet 
quenching in final state interaction have to be considered. 
Jet quenching in gen eral suppresses high transverse mo- 
mentum hadrons [16[ . If we assume the effects of parton 
and hadron final state interactions on the total hadron 
multiplicity to be negligible fl7l - [i"9| , the only uncertainty 
for hadron multiplicity density in A + A collisions comes 
from the nuclear modification of parton distribution func- 
tions at small x. HIJING2.0 employes the factorized form 
of parton distributions in nuclei, 



tf(x,Q')=ARA(x,Q')f«(x,Q') 



(1) 



where R^(x,Q 2 ) is nuclear modification factor as given 
by the new HIJING parameterization Q, 
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FIG. 1: (color online) HIJING2.0 parameterization of nuclear 
shadowing factor for light quark distributions in lead nuclei 
(solid) with s q = 0.1 as compared to EPS09 (dashed and 
dotted) [U and HKN07 (dot-dashed and long dashed) [HI 
parameterizations at two different scales. 
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Rf(x,b) = 1.0 + 1.19 log 1/( k (a; 3 - 1.2.T 2 + 0.21x) 

- s q {b) {A 1/3 - 1)°- 6 (1 - 3.5V5) 

x exp(-x 2 /0.01), (2) 

Rf(x,b) = 1.0 + 1.19 log 1/( k (x 3 - 1.2.t 2 + 0.21x) 

- s g (b) (A 1 / 3 - lf 6 (l - 1.5x ' 35 ) 

x exp(-a; 2 /0.004), (3) 

for quarks and gluons, respectively. The impact- 
parameter dependence of the shadowing is implemented 
through the parameters, 



FIG. 2: (color online) The same as Fig. Q] except for gluon 



s a (b) = Sa -(l-b 2 /R 2 A ), 



(4) 



where Ra = 1.1 2 A 1 / 3 is the nuclear size. The form of 
the impact parameter dependence is chosen to give rise 
to the centrality dependence of the pseudorapidity multi- 
plicity density per participant pair 2dN ~ c h / 1 dr\ / '-/Vp art . The 
impact-parameter averaging is done with a weight of the 
thickness function of the nucleus with a hard-sphere nu- 
clear distribution. 



distributions with s g = 0.22 — 0.23. 



The assumed form of nuclear modification in Eq. (fT]) 
has been studied with data from deeply inelastic scatter- 
ing (PIS) and Drell-Yan lepton pair production experi- 
ments [lo, Hi]]. I n the new HUNG parametrization, the 
value s q = 0.1 is fixed by the experimental data on DIS 
off nuclear targets 0. Shown in Fig. Q] is the HIJING2.0 
nuclear shadowing factor for light quarks {u and d) as 
compared to that by EPS09 and HKN07 [U pa- 
rameterizations at two different scales Q — 2.0 and 4.3 
GeV/c. Both EPS09 and HKN07 use vacuum DGLAP 
evolution equations for the parton distributions to de- 
termine the scale dependence of the nuclear shadowing 
factors which is rather weak for the quark distributions. 

The gluon shadowing at small x is not constrained di- 
rectly by DIS and Drell-Yan experimental data, except 
momentum conservation. Total momentum depletion 
due to suppression of both quark and gluon distributions 
at small x is partially compensated by the anti-shadowing 
at x w 0.1 for quarks. Momentum conservation, however, 
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will not provide much constraints on gluon distribution 
in nuclei both in the shadowing and anti-shadowing re- 
gion as shown by EPS09 and HKN07 paramctcrizations 
in Fig. [3J both enforce momentum conservation. HI- 
JING2.0 parameterization violate momentum sum rule 
by about 16% which can be easily compensated by ad- 
justment of gluon modification factor at large and inter- 
mediate x which will not affect mini-jet production at 
small x. 

The value of gluon shadowing parameter s g in HI- 
JING2.0 is constrained only by the hadron multiplicity 
in heavy-ion collisions. Using the combined RHIC data 
[H on the centrality dependence of charged hadron mul- 
tiplicity density in mid-rapidity as constraints, a range 
s g = 0.17 - 0.22 was obtained 0]. Note that HI- 
JING2.0 parameterization assumes a scale-independent 
form. Such an approximation for gluon distribution will 
not be valid at very large scale due to dominance of gluon 
emission dictated by the DGLAP evolution equations. 
Both EPS09 and HKN07 paramctcrizations use the vac- 
uum DGLAP evolution equations for parton distribu- 
tions to determine the scale dependence of the nuclear 
shadowing factors. At Q = 2.0 and 4.3 GeV/c, which 
are the typical scales for minijet production at RHIC 
and LHC, respectively, the gluon shadowing varies at 
most about 13% in EPS09 parameterization as shown in 
Fig. [2J As shown in the figure, this is also approximately 
within the uncertainty in the HIJING2.0 parameteriza- 
tion of gluon shadowing constrained by the charged mul- 
tiplicity data in heavy-ion collisions (s g = 0.22 — 0.23). 
Furthermore, higher-twist contributions to the DGLAP 
evolution equations might become important and would 
modify the scale dependence of the shadowing factor at 
intermediate scales. 

With the above parameterization of parton shadow- 
ing and the range of gluon shadowing parameter s g = 
0.17 — 0.22, the predicted 2dN c h/d?]/N pm - t , shown as 
dash-shaded area in Fig. [31 agrees well with the new 
ALICE data in the most central Pb + Pb collisions at 
y/s = 2.76 TeV, within the experimental error and a 
large theoretical uncertainty of about 15% from the gluon 
shadowing parameter. The HIJING2.0 results are ob- 
tained by calculating dN c h/dr] and V part for different 
impact-parameters squared b 2 with equal bin size. By 
performing a combined x 2 -fit of the RHIC data [l| for 
Au + Au collisions at ^fs = 200 GeV and the data point 
from ALICE in the most central Pb + Pb collisions at 
\/s = 2.76 TcV, the range of gluon shadowing parame- 
ter is reduced to s g = 0.20 — 0.23. As shown in Figs. Q] 
and [21 the HIJING2.0 gluon shadowing from such fit is 
much stronger than the EPS09 and HKN07 parameteri- 
zations. However, it is comparable to the parameteriza- 
tion by Strikman et al. [22[ and both are much stronger 
than the nuclear shadowing for quark distributions. 

With this new range of s g and therefore reduced un- 
certainty we calculate the prediction for the centrality 
dependence of dN c h/drj in Pb + Pb collisions at both 
-^/s = 2.76 and 5.5 TcV, shown in Fig. [31 as solid-shaded 
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FIG. 3: (color online) Charged hadron multiplicity density 
in mid-rapidity per participant pair 2dN c h/dr]/Np al t as a 
function of iV pa rt from HIJING2.0 calculation with gluon 
shadowing parameter s g = 0.20 — 0.23 (solid-shade) and 
s g — 0.17 — 0.22 (dash-shade) as compared to combined RHIC 
data [lj for Au + Au collisions (filled circle and star) and AL- 
ICE data [j| at LHC (solid square). 
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FIG. 4: (color online) The ratio of charged hadron multiplicity 
density in mid-rapidity in heavy-ion collisions at different col- 
liding energies, using Au + Au collisions at yfs = 200 GeV as 
the common denominator. The data at the RHIC energies are 
combined from different experiments [lj. The ALICE data on 
non-single diffractive (NSD) pp collisions at yfs = 2.36 TEV 
[23l ] is used to get the data point at y/s — 2.76 TeV using HI- 
JING2.0 extrapolation. The NSD pp data at s/s = 0.2 TeV 
is from UA1 Hi 



area. The calculated centrality dependence of dN c h/dr] 
in Au + Au collisions at two RHIC energies is also shown 
together with combined RHIC data [l|. 

To exam the centrality dependence of dN c h / drj at dif- 
ferent colliding energies in detail, we plot in Fig. 2] the 
ratio of dN c h/drj at different colliding energies using 
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Au + Au collisions at s/s = 0.2 TeV as the common de- 
nominator. In the figure we also plot (pp 2.76 TeV)/(pp 
0.2 TeV), using the ALICE data on non-single diffractive 
pp collisions at y/s = 2.36 TeV [H and HIJING2.0 cal- 
culation to extrapolate to the value at y/s = 2.76 TeV. 
The pp data at ^fs = 0.2 TeV is from UAI [H. The 
ratios of charged hadron multiplicity densities at the two 
LHC energies to that at RHIC are surprisingly flat over 
a large range of centralities just as the ratio of two RHIC 
energies. It is interesting to note that the increased en- 
ergy dependence of charged multiplicity density in central 
Pb + Pb collisions over that in pp is reached already at 
iV part = 50(100) for sfs = 2.76 (5.5) TeV. In other words, 
the centrality dependence of charged hadron multiplicity 
density increases with energy in peripheral Pb + Pb col- 
lisions. Such centrality dependence of charged hadron is 
a consequence of the impact-parametcr-dependent gluon 
shadowing in HIJING2.0. 

With a given transverse momentum cut-off pg, the to- 
tal number of mini-jets per unit transverse area could 
become so large that it exceeds the limit 

T AA (b)a jet pi 

ttR\ ~ ~ (5) 

for independent multiple jet production for sufficiently 
large inclusive jet cross section at high colliding ener- 
gies and for large nuclei, where T A a(V) is the overlap 
function of A + A collisions and tt/pq is the intrinsic 
transverse size of a mini-jet with transverse momentum 
Pq. This is the reason for an energy-dependent cut-off 
Po for high-energy pp collisions in HIJING2.0 since the 
GRV parton distributions [r| have a large gluon dis- 
tribution at small x and therefore large mini-jet cross 
section at high colliding energies. The above limit for in- 
coherent mini-jet production should also depend on nu- 
clear size and impact-parameter which can be determined 
sclf-consistently through Eq. © QJ]. In HIJING2.0 such 
impact-parameter dependence of the cut-off scale is not 
considered. Instead, an impact-parameter dependence 
of the gluon shadowing in Eq. ^ is considered that is 



stronger than the typical nuclear length La = \J R\ — b 2 
dependence. Such a stronger impact-parameter depen- 
dence is favored by the centrality dependence of dN c h/r] 
in Au + Au collisions at RHIC. This is also the reason for 
nearly energy-independence of the centrality dependence 
of charged hadron multiplicity density at the LHC ener- 
gies. If confirmed by experimental measurements, it will 
have important implications on the initial eccentricity for 
the study of elliptic flow and jet quenching. 

In summary, we have carried out a combined fit of the 
new ALICE data Q on charged hadron multiplicity den- 
sity in the most central Pb + Pb collisions at ^/s = 2.76 
TeV and the RHIC data within HIJING2.0 model. The 
range of gluon shadowing parameter s g = 0.20 — 0.23 
in the new HIJING parameterization of parton shadow- 
ing [1] enables us to predict the centrality dependence of 
the charged hadron rapidity density with reduced uncer- 
tainty in Pb+Pb collisions at y/s = 2.76 and 5.5 TeV. The 
centrality dependence is surprisingly independent of col- 
liding energy for most centralities starting at Apart = 50 
(100) for -^=2.76 (5.5) TeV. However, the centrality de- 
pendence in the peripheral collisions becomes stronger at 
higher colliding energies. 
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